How Does Casimir Energy Fall? 
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Doubt continues to linger over the reality of quantum vacuum energy. There is some question 
whether fluctuating fields gravitate at all, or do so anomalously. Here we show that for the simple 
case of parallel conducting plates, the associated Casimir energy gravitates just as required by the 
equivalence principle, and that therefore the inertial and gravitational masses of a system possessing 
Casimir energy E c are both E c /c 2 . This simple result disproves recent claims in the literature. We 
clarify some pitfalls in the calculation that can lead to spurious dependences on coordinate system. 
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The subject of quantum vacuum energy (the Casimir 
effect) dates from the same year as the discovery of renor- 
malized quantum electrodynamics, 1948 It puts the 
lie to the naive presumption that zero-point energy is 
not observable. On the other hand, it continues to be 
surrounded by controversy, in large part because sharp 
boundaries give rise to divergences in the local energy 
density near the surface (see Refs. @, H, The most 
troubling aspect of these divergences is in the coupling 
to gravity. Gravity has its source in the local energy- 
momentum tensor, and such surface divergences promise 
serious difficulties. The gravitational implications of 
zero-point energy are an outstanding problem in view of 
our inability to understand the origin of the cosmological 
constant or dark energy [5|, [6|, [7| . 

As a prolegomenon to studying such issues, we here ad- 
dress a simpler question: How does the completely finite 
Casimir energy of a pair of parallel conducting plates cou- 
ple to gravity? The question turns out to be surprisingly 
less straightforward than one might suspect! Previous 
authors 

[I H M, ED, El have given disparate answers, 
including gravitational forces, or gravitationally modi- 
fied Casimir forces, that depend on the orientation of the 
Casimir apparatus with respect to the gravitational field 
of the earth. We will here resolve some of this confusion 
with a convincingly calculated result consistent with the 
equivalence principle. That is, the renormalized Casimir 
energy couples to gravity just like any other energy. In 
our opinion, this fact is evidence that vacuum energy 
must be taken seriously in gravitational theory and that 
the problem of boundary divergences must be resolved by 
a better understanding of the modeling and renormaliza- 
tion processes. 

We start by recalling the electromagnetic Casimir 
stress tensor between a pair of parallel perfectly con- 
ducting plates separated by a distance a, with transverse 
dimensions L»a, as given by Brown and Maclay 



where the third spatial direction is the direction normal 
to the plates. This is given in terms of the Casimir en- 
ergy per unit area, £ c = — ir 2 he /(720a 3 ). Outside the 
plates, (T 1 *") = 0. Omitted here is a constant divergent 
term that is present both between and outside the plates, 
and also in the absence of plates, which cannot have 
any physical significance. Because the electromagnetic 
field respects conformal symmetry, there is no surface di- 
vergent term such as is present for a minimally coupled 
scalar field subject to Dirichlet conditions on the plates, 
or more generally for curved surfaces [3]. (Henceforth 
we will set ft = c = 1.) 

Now we turn to the question of the gravitational in- 
teraction of this Casimir apparatus. It seems to us that 
this question can be most simply addressed through use 
of the gravitational definition of the energy-momentum 
tensor, as a variation of the matter part of the action, 



W a = SW„ 



(2) 



Following Schwinger (note the factor of 2 in the defini- 
tion), for a weak field we define g^ v — ri^ v -\-2h tlll . To first 
order we can ignore y/—g. The gravitational energy, for 
a static situation, is therefore given by (SW = — f dt 8E) 



E g = SE m = - / (dx) h^T^ 



(3) 



We then replace here by the one-loop expectation 
value of the electromagnetic stress tensor ([T]). 

Calloni et al. 0] and Bimonte et al. [13] use the Fermi 
metric 



5oo 



(1 + 2gz) , gij = 6, 



(4) 



(7**) =^diag(l,-l,-l,3), 
a 



(1) 



in terms of the gravitational acceleration g. This is evi- 
dently appropriate for a constant gravitational field. We 
will discuss its relation to the field due to the earth below. 
Let the apparatus be oriented at an angle a with respect 
to the direction of the gravitational field. The Cartesian 
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FIG. 1: Relation between two Cartesian coordinate frames: 
One attached to the earth (x,y, z), where — z is the direction 
of gravity, and one attached to the parallel-plate Casimir ap- 
paratus (C, fj, x)> where £ is in the direction normal to the 
plates. The parallel plates are indicated by the heavy lines 
parallel to the r\ axis. The x — \ axis is perpendicular to the 
page. 



coordinate system attached to the earth is denoted by 
(x, y, z), where, as noted above, z is the direction of — g. 
The Cartesian coordinates associated with the Casimir 
apparatus are (Ci^jX): where C is normal to the plates, 
and rj and \ are parallel to the plates. (See Fig. [TJ) The 
center of the apparatus is located at (Co , 77 = 0, \ — 0). 
Now we calculate the gravitational energy from Eq. ((3j) : 

a 



E g = J {dx) gzT 00 = — L 2 a Co cos a, 



(5) 



up to an additive constant, independent of Co ■ (Any 
constant added to the gravitational potential gz simply 
shifts the value of this constant.) Thus, the gravitational 
force per area A = L 2 on the apparatus is 



-Be 



(6) 



F _ 1 dE g 

~A ^~A ~dz^ 3 ~ u ~ 2a 

because zq = Co cos a. Note that on the earth's surface, 
the dimcnsionlcss number e = 2ga/c 2 is very small. For 
a plate separation of 1/xm, e = 2 x 10~ 22 , so the effects 
that we are considering are very small compared to the 
gravitational forces on the plates. (However, Calloni et 
al. discuss the possibility of experimentally observing 
this force.) 

It is a bit simpler to use the energy formula ((3]) to 
calculate the force by considering the variation in the 
gravitational energy directly, that is, 



5E g = - I (dx) Sh^T^. 



(7) 



It is easy to verify that this gives the correct force on a 
mass point, F = mg. If we use this formula to calculate 
the gravitational force on the rigid Casimir apparatus, by 
considering a virtual displacement upward by an amount 
8zq, we find the same a-independent result (|6|). 

Alternatively, we can start from the definition of the 
gravitational field [l5| , 



SW g = 



J(dx) 



ST^h 



(8) 



which can again be checked to yield the correct force on 
a mass point. For the constant field ([4} the force on a 
Casimir apparatus is obtained from the change in the 
energy density T 00 ; that is, recalling that zq = Co cos a, 



_ £cSz a 1 



-[«5(C-Co-a/2)-<5(C-Co+a/2)], (9) 



a cos a 

where the S functions arise from the step functions at the 
boundaries. This yields from Eq. (JSj) the same result J6J. 

Our answer is consistent with the principle of equiva- 
lence, and with the second analysis of Jaekel and Rey- 
naud fl6| . who state that the inertia of Casimir energy 
(at least in two dimensions) is E c /c 2 . However, it is only 
i that found by Bimonte et al. [l2:], which is also the 
first force formula [Eqs. (7) and (8)] provided by Cal- 
loni et al. Q. Our Eq. does, however, reproduce the 
second formula [Eq. (9)] given in Ref. [9j, which those 
authors describe as the one that should be observable. 
We discuss this situation further below. 

We now digress to consider whether the constant- 
field approximation ^ is adequate for an apparatus 
suspended above the earth or some other pointlike 
mass. Should we instead use the perturbation of the 
Schwarzschild metric? One might expect that the result- 
ing curvature corrections are of order ^ < 1, at worst, 
relative to the main term, where R is the earth's radius. 
The point, however, is that naive attempts to do the cal- 
culation in curved space change the answer by factors like 
2, and also differ among themselves, and the resolutions 
of the fallacies are sufficiently instructive to justify our 
belaboring the point. 

The Schwarzschild metric in isotropic coordinates 17 1 
and for weak fields (GM/r <C 1) is 



/ 2GM\ , 2 
- - ' 1 — ) dt 2 



ds 2 = - l- 



1 



2GM\ 



dv 2 



(10) 



If we expand a short distance z above the earth's sur- 
face, we find the nonzero components of the gravitational 
field to be hoo = hu = I122 = /133 = GM/R — gz, in 
terms of the acceleration of gravity, g = GM / R 2 . (It 
is important to recognize that the constant GM/R is ir- 
relevant in the following, and that correspondingly the 
results do not depend on the origin of z.) The virtue 
of isotropic coordinates is that the spatial line element 
(apart from an overall factor) has the usual Cartesian 
form dr 2 = dx 2 + dy 2 + dz 2 . Now when we compute the 
gravitational energy from Eq. ([7|) each component of the 
Casimir stress tensor contributes with equal weight: 

Ka+a/2 

SE g = gASz ddT^ + ^+T^ + T 33 ), (11) 

J Co -a/2 

which gives the force 
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-2gS c 



(12) 
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since T = T x \ = 0. This is twice the previous result ([6]). 
Note that again the result is independent of a. The same 
result is obtained if we start from Eq. (jHJ). 

We should be able to obtain the same result using 
the original Schwarzchild coordinates, where ft-oo = — 9 Z -, 
h pp = —gz, and all other components of are zero. 
However, now if we use the first method ([7]), the result 
is F/A = —4g£ c cos 2 a, so now the force depends on the 
orientation of the apparatus. Even if a — 0, the magni- 
tude differs from Eq. (I12|) by an additional factor of 2. 
(It thereby fortuitously agrees with the result in Ref. [l2| 
for that angle.) 



What is going on here? The reason we get differ- 
ent answers in different coordinate systems is that our 
starting point (|3|) is not gauge-invariant. Under a co- 
ordinate redefinition, which for weak fields is a gauge 
transformation of [15|, — > h^ v + d^ v + , 
where £p is a vector field, Eq. ([3]) is invariant only if 
the stress tensor is conserved, d^T^ — (in the weak- 
field context). Otherwise, there is a change in the action, 
AW g — —2 j(dx)£ :l/ d tl T tiU . Now in our case (where we 
make the finite size of the plates explicit, but ignore edge 
effects on T^ v because L a) 



T K» = ^ diag(l, -1, -1, 3)0(C - Co + a/2)6(a/2 - Q + ( a )9( V + L/2)6(L/2 - 7y)0( X + L/2)6(L/2 - X ). (13) 



Taking the divergence of Eq. (|13jl gives corresponding S functions on the surfaces and leads immediately to 

6£ f IE r 

AE g = -^j dr 1 dxMCo-a/2, V , X )-^(Co + a/2, V , X )]~^ J d^d X [£„(£ -L/2, X ) ~ £„(C, L/2, X )] -(»?«-> X}- 



This transformation entirely accounts algebraically for 
the difference between the forces in isotropic and 
Schwarzchild coordinates, but it does not yet explain 
physically why there arc two different answers, nor tell 
us which, if either, is correct. 

There seem to be two possible ways to proceed. First, 
it is clear that the energy-momentum tensor of the com- 
plete physical situation must be conserved, and there- 
fore the expression ([3]) would be gauge-invariant if we 
included a physical mechanism holding the plates apart 
against the Casimir attraction. That road probably leads 
to complicated, model-dependent calculations. The al- 
ternative is to find a physical basis for believing that one 
coordinate system is more realistic than another. Fortu- 
nately, that problem apparently has a natural solution. 
The crux of the difficulty is that the relations between co- 
ordinate increments and physical distances depend upon 
the distance from the gravitating center in the most com- 
mon coordinate systems. 

Of course, a perfect coordinate system is not possible 
in a curved space, but the kind that comes closest to rep- 
resenting distances accurately all along a timelike world- 
line is a Fermi coordinate system, the general-relativistic 
extrapolation of an inertial coordinate frame. Such a 
system has been given by Marzlin [l8| for a resting ob- 
server in the field of any static mass distribution. Here we 
give a simple rederivation for the case at hand. Starting 
from the isotropic metric (|10|) , first eliminate the con- 
stant term by rescaling the coordinates, t — > (l + ^jp) t, 



r — * (l — ^jjp) r, and expand to first order: 

ds 2 = -(1 + 2gz)dt 2 + (1 - 2gz)dr 2 . (15) 

But we need r to measure physical displacements even 
when z ^ 0, so we write 



x = x' + gx'z', y = y' + gy' 

2 1 



z = z' + Uz' 2 



x' 2 -y' 2 ) 



(16) 



Then to first order in coordinates we obtain the Fermi 
metric ((4]). The corresponding gravitational force is 
therefore given by Eq. ([6]), after all! 

Now we can use the method described above to trans- 
form the energy in isotropic coordinates to that in Fermi 
coordinates. We use Eq. (fT4|) to compute the additional 
gravitational energy in terms of the gauge field £ M that 
carries us from isotropic coordinates to Fermi coordi- 
nates, 



(17) 



Here h 1 ^ = -gz, h\ } = -gzSij , h£ = -gz, hf 3 = 0, 



IF 

h Q l = 0. The gauge field turns out to be 



1 / 1 



-5 -( cos a + (rj sin a 



1 



1 



—g £77 cos a H — 77 sin a 



f(ri,X), 

g((,x), 



(Ccosa + 77 sin a) X + h{C,,rf), 



(18) 
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where the functions /, g, and h are irrelevant. Now from 
Eq. (TT4]) we obtain a AE g (zo) that yields an additional 
force, AF/A = g£ c . When this is added to the isotropic 
force (fT2"]) , we obtain the result given in Eq. © 

F 1 + AF F F 

^ = -2g£ c + g£ c = -g£ c = — . (19) 

The conceptual reason why other coordinates give dif- 
ferent answers is that under the virtual displacement in- 
volved in defining one is stretching the apparatus as 
well as moving it. Correcting for the spurious changes in 
L and a restores the Fermi result in all cases. The im- 
portance of distinguishing a from the physical gap was 
noted by Sorge llj in studying a related problem. 

As noted above, Calloni et al. [9J find a result 4 times 
ours, which is the only result from Ref. cited in the 
later paper I2| with which it shares authors. How- 
ever, Ref. [9] states that that force formula has two 
parts, in the ratio of 3/1, and that only the smaller 
piece is "Newtonian," or "to be tested against observa- 
tion." Our understanding of what that means is the 
following. Start from Eq. ((2]) and consider a general 



coordinate transformation, x'^ = 
9uu( x ') = 9^{x) + Sg^(x), where 



x^ 1 + Sx^, so that 



5g^ y = 5x x d\g t _ LU + < 



d6x a 
' dx^ 



AW? 



dSx? 
dx v 



(20) 



For a rigid translation, Sx x is a constant, so only the first 
term in Eq. (|20|) is present, which gives the result |6]). 
However, if we do not make this restriction, we obtain 
from Eq. ^ (after integration by parts) a surface-term 
correction to the force: 



-9 fx = 



SW m 

Sx x 



on 



(21) 



where the force vector density is 0] —V„T U \ or 



9T v x)+-sf=-gT» v d x9tlv . (22) 



Note that the surface term identically cancels the first 
term in f\ . Now if Q, refers to all space, the surface 
term vanishes (as we have shown explicitly). But if 
is just the space volume between the plates, and we in- 
clude this correction for the Fermi metric ^ for which 
\J—g = 1 + gz, we obtain an additional term — 3g£ c cosa. 
Adding this to the previous result ©, we obtain the re- 
sult of Ref. [13] if a = 0. However, in general the result 
depends on the angle between the apparatus and the ver- 
tical. Is this consistent with the equivalence principle 
(the scalar nature of mass)? A similar angle dependence 
will now occur with the isotropic Schwarzchild metric. 
Why should one trust the formula (|22p over the more 
fundamental variational principle when boundaries are 
present? Omission of the surface term resolves the dis- 
crepancy, giving the equivalence-principle result (0). 



After circulation of the present paper as 
|arXiv:hep-th/070209l| Bimonte et al. responded 
that (i) our result is implicit in their earlier analysis; 
(ii) our Eq. @ implicitly assumes the equivalence 
principle, which they identify with the conservation law 
V^F" = 0. We reply: (i) Ref. [H states without quali- 
fication that the force is 4 times the correct value, (ii) 
By "equivalence principle" we merely mean that gravity 
tends to enforce geodesic motion. This is not obviously 
equivalent to conservation, as shown by the existing 
uncertainty in the literature. Contrary to statements 
in Ref. [19(, it is generally accepted in quantum field 
theory in curved space-time that the renormalized total 
T^ v must be conserved. 
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